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THE DISSIPATION RATE TRANSPORT EQUATION 
AND SUBGRID-SCALE MODELS IN ROTATING TURBULENCE* 

ROBERT RUBINSTEIN^ AND YE ZHOU* 

Abstract. The dissipation rate transport equation remains the most uncertain part of turbulence mod- 
eling. The difficulties arc increased when external agencies like rotation prevent straightforward dimensional 
analysis from determining the correct form of the modelled equation. In this work, the dissipation rate 
transport equation and subgrid scale models for rotating turbulence are derived from an analytical statisti- 
cal theory of rotating turbulence. In the strong rotation limit, the theory predicts a turbulent steady state 
in which the inertial range energy spectrum scales as k~ 2 and the turbulent time scale is the inverse rotation 
rate. This scaling has been derived previously by heuristic arguments. 

Key words. Dissipation rate transport equation, subgrid scale models, rotating turbulence 

Subject classification. Fluid Mechanics 

1. Introduction. Kraichnan’s (1959) Direct Interaction Approximation (DIA) and related Lagrangian 
closures (Kraichnan, 1964; Kaneda, 1968) remain the only fully deductive turbulence theories. Although a 
study of a simple inhomogeneous flow, like channel flow, using these closures would be of the greatest theo- 
retical and practical interest, the complexity of the calculations required has precluded any but preliminary 
results (Danncvik, 1992). 

Practical application of DIA therefore requires some compromise of rigor in the interest of utility. The 
most comprehensive attempt to extract turbulence models from DIA remains the two-scale theory (TSDIA) 
of Yoshizawa (1984, 1996) in which inhomogeneity, anisotropy, and time- dependent nonequilibrium effects 
are introduced by perturbing about a state of homogeneous, isotropic, stationary turbulence. 

Yoshizawa’s procedure leads to formulas for quantities familiar in single point phenomenological turbu- 
lence closures like the two-equation model. A typical result (Yoshizawa, 1984) is the expression for eddy 
viscosity 

(i) ^ = ^ j dk r dr r) 

in terms of the DIA descriptors (Kraichnan, 1959) of isotropic turbulence: the response function G(k ) r) 
and correlation function Q(/c,r), The isotropy of the lowest order field implies that these descriptors are 
scalars, homogeneity permits introduction of the wave-vector argument k, and stationarity in time makes 
them functions of time difference r only. Yoshizawa (1984) shows how the familiar eddy viscosity formula 
of single-point turbulence modeling is deduced from this formula, by substituting Kolmogorov scaling forms 
for G and Q . As is well-known, Kolmogorov scaling describes the simplest turbulent steady state with a 
constant flux of kinetic energy from the large to the small scales. 

2. DIA for Rotating Turbulence. Eq. (1) suggests that an eddy viscosity for turbulence subject to 
any external agency can be derived, provided that appropriate formulas for G and Q are known. Rotation 
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is a particularly simple external effect, since the energy remains an inviscid invariant under rotation, and a 
steady state with constant energy flux remains possible. But Kolmogorov scaling may no longer apply to 
this steady state and the occurence of a distinguished time scale precludes the deduction of the applicable 
scaling law by dimensional analysis. To deduce the scaling, we will appeal to closure in the form of the direct 
interaction approximation. 

For rotating turbulence, the DIA equations of motion take the form 


( 2 ) 


(3) 


Gij (k, t, s) + t 2Pip(k)£lpqGqj (k, t , s) 
+ J dr T)ip(k,t,r)G P j(k,r,s) =0 
Qij (k, £, s) -h 2 P ip (k)f7p <? (5gj (k, t, s) 

+ Is dr r? ’ p ( k ’*’ r ^ p: '( k ’ r ’ S ) 

— f dr G jp(k, t , v^Fpj (k, s, 7™) 

Jo 


where the eddy damping 7] and forcing F are defined by 


(4) 

( 5 ) 


77ir(k, t)S} — / ^pdc|Fj mn (k)F jUrs (p) x 

J k=p+q 

Gmp, (p, ty s)Qns(Q) 5 ) 

Fij(k,t,s) = f d , pdciPi Tnn (]&)Pj rs (k) x 
7k— p+q 

Qns (P? S)Qmr (Q> f j 


In Eqs. (4) and (5), 

Tjrnn(h) = kmPinO*) + k n P xrn { k) 

P%m (k) = di m ~ k 

and ftpq is the antisymmetric rotation matrix. The solution of these equations in complete generality is 
not known. A useful simplification, EDQNM, effectively replaces the response equation Eq. (2) by a 
phenomenological hypothesis, and solves a simplified, Markovianizcd version of the correlation equation Eq. 
(3) (Cambon and Jacquin, 1986). 

A perturbative solution of these equations is suggested by Leslie’s (1972) treatment of shear turbulence, 
treat the rotation terms as small, and perturb about an isotropic turbulent state. This approach is adopted 
by Shimomura and Yoshizawa (1986) who derive a TSDIA theory in which both inhomogeneity and rotation 
are described by small parameters. 

A complementary limit is also of interest. Namely, in the response equation, balance the time derivative 
by the rotation term, and treat the eddy damping as small. This linear theory of the response equation 
treats strongly rotating turbulence as a case of weak turbulence (Zakharov et al , 1992) in which nonlinear 
decorrelation of Fourier modes is dominated by linear dispersive decorrelation (Waleffe, 1993). The result is 
conveniently expressed in terms of the Craya-Herring basis 


e^(k) — kxrt/ | kxfl | 
e (2) (k) = kx(k xfi)/ | kx(kxft) | 
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or the equivalent basis of Cambon and Jacquin (1989), and the corresponding tensors 


,0 =e (l) c (2)_ (1) (2) 

^ = ef)ef + e«ef> 

f2 _ M)J1) 


( 6 ) 


p. = e S l> e K 1 


£ 3 . = e- 


(2)^(2) 


e) + e) 


Note that £?- = Pij(k). 

The leading order solution of Eq. (2), obtained by dropping the eddy damping term, is 


Gij( k, t , s) = {cos(k ■ rt(t — s)/k)Pij( k) 

(7) + sin(k Sl(t- $)/k)§ (k )}H(t - s) 

where H is the unit step function. Let us adopt the fluctuation- dissipation hypothesis relating the two-time 
correlation function to the response function and single- time correlation function 

(®) Qij (k> ^ s ) = G\m{ k 5 s)Q rn j (k) + (j jm (k, s, i)Qmi{ k) 

Conditions under which this approximation is reasonable are discussed by Woodruff (1992). Substituting 
Eqs. (7) and (8) in Eq. (3) shows that the single-time correlation function must take the general form 
containing all of the £ tensors of Eq. (6) 

(®) Q.,(k) = ^Q p (k)^(k) 

which is equivalent to the form of the correlation function noted by Cambon and Jacquin (1989). 

The DIA inertial range energy balance (Kraichnan, 1971), which states that a steady state with constant 
energy flux exists, is 

r OC 

e=[I + - I~]Pimn (k) / dr 2P firs (p)G rnfl (p y r)Q ns {q, r)Q ir (k, r) 

Jo 

(1®) — PjrsO*)Gij (k, r)Q ns ( p, r)Q mr (q, r) 

where the integration operators in Eq. (10) are defined by 


/ + (*o) = 

/ dk 

1 dpdq 


J k>k 0 

/k=p-f-q;p 1 q<&o 

I~(k 0 ) = 

f dk 

1 dpdq 


Jk<k 0 

/k=p+q \p,q>k 0 


and r = t~s denotes time difference. The time integrals in Eq. (10) will take the form 

(ii) e(k, p, q ) = n~ l s(± P ■ n/ P n ± q . n/ q n ± k • n/*n) 

where 


SI = (S2 pg n pg y/ 2 

Thus, wave-vector integrations take place over resonant triads only (Waleffe, 1993) and these integrals scale 
as Sl~ l (Zhou, 1995). 
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Introduce the ansatz 


(12) Q p (k) = k~ a ~ 2 f p (k, n ■ k/fcfi) 

for the functions in Eq. ( 9 ). In Eq. ( 12 ), f p is homogeneous of degree zero in k. In view of Eq. ( 11 ), the 
energy balance Eq. ( 10 ) requires a = 2, as in the heuristic argument of Zhou (1995). The goal of this work 
is a two-equation model of the standard form. Accordingly, we only attempt to quantify the overall effect of 
rotation on turbulent energy transfer, ignoring the polarization of turbulence (Cambon and Jacquin, 1989) 
by rotation and the distribution of energy in k space. For this purpose, it suffices to replace the anisotropic 
energy spectrum tensor by the energy transfer-equivalent isotropic spectrum given by Zhou (1995), 

(13) E(k) ~ \Zeflk~ 2 


It is of interest to exhibit the eddy damping correction to the leading order results represented by Eqs. 
( 7 ) and (13). Namely, substituting these results in Eq. (4) for the eddy damping factor, 

roc ^ 

(14) 77 ~ /c 2 / dp —\feVtp ~ 2 ~ ky/efVt 

Jk * * 

The corrections to the time scale and energy spectrum have the form 


(15) 


e ~ ^{ 1 + °( n-3/2 )}” 1 

E ~ v'en*T 2 { 1 + 0((fc 2 e/0 3 ) 1/2 } 


The low rotation rate expansion of Shimomura and Yoshizawa (1986) gives the complementary expansions 
in positive powers of fi, 

e~e- 1 ' 3 k- 2 t 3 {l + 0(Cl)} 

(16) E ~ e 2/3 k- 5/3 {l + 0(tl/e 1/3 k 2/3 } 

The expansions in Eqs. (15) and (16) could be consolidated into Padc approximations for the energy spectrum 
and decorrelation time applicable for intermediate rotation rates. 

2 . 1 . Locality of energy transfer. The scaling law of Eq. (13) is purely formal: to prove that a 
steady state scales this way, the convergence of the integral in Eq. (10) must be demonstrated. Divergence 
would imply strong dependence on the cutoff at large or small scales, and would therefore alter the scaling 
law (Kraichnan, 1959). Even if the dependence of E(k) on ft 1 / 2 is known, Eq. (13) cannot be asserted on 
dimensional grounds: dimensional analysis assumes that k~ l is the only relevant length scale; it therefore 
postulates the locality which convergence of the flux integral proves. 

However, it is not difficult to prove that the flux integral does in fact converge when evaluated for an 
infinite Ar ~ 2 inertial range. The proof only requires the transverse character of the tensors £ p and does not 
differ conceptually from the proof originally given by Kraichnan (1959). Accordingly, we can assert the 
locality of energy transfer in a rotating inertial range and can apply the scaling laws of Eq. (15) to develop 
turbulence models. 

3 . Subgrid-Scale Models. Since the lowest order field in TSDIA is arbitrary, the transport properties 
of weakly strained strongly rotating turbulence can be derived by perturbing about the steady rotating 
turbulent state just described instead of perturbing about a Kolmogorov steady state. In carrying out this 
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program, it should be stressed that formulas like Eq. (1) for the eddy viscosity are the result of evaluating 
an integral over interacting triads which satisfy a resonance condition. This condition provides an additional 
source of anisotropy in the exact theory of rotating turbulence, since only two-dimensional triads for which 

k n = p fi = q n = o 

satisfy the resonance condition automatically. However, since our goal is to evaluate the over-all effect of 
rotation on turbulent energetics, we will evaluate the eddy viscosity integral using the isotropic expressions 
ofEq. (15). 

The simplest way to derive a sub-grid scale viscosity from Eq. (1) is to integrate only over scales 
satisfying k > 2n/A, where A denotes the filter width. It follows from substitution of Eq. (15) in Eq. (1) 
that in strongly rotating turbulence, to leading order 

( 17 ) v = C?A^/n 

where A is the filter size. Corrections for finite rotation rates are suggested by Eq. (15). Setting A equal to 
the integral scale of turbulence leads to the single point result 

(18) 

n 

A subgrid scale model is derived by eliminating the dissipation rate from Eq. (17). We follow the 
derivation of the Smagorinsky model by equating the dissipation rate to the resolved production: 

( 19 ) e= l -uS lj S ij 
where S 2 = S XJ Sij/2 and Sij is the resolved strain rate. Therefore, 

(20) £ = C?y/eJtiAS 2 
Solving Eq. (20) for e and substituting in Eq. (17), 

(21) „ . <C?) J ^ 

which is the Smagorinsky model for strongly rotating turbulence. 

This calculation supplements the evaluation by Shimomura and Yoshizawa (1986) of the sub-grid scale 
viscosity in weakly rotating turbulence. But the effect of rotation becomes weaker with decreasing scale size; 
accordingly, the weak rotation correction of Shimomura and Yoshizawa may be appropriate even in some 
rapidly rotating flows. 

Additional production and dissipation mechanisms are often modeled by Richardson number modifica- 
tions of the subgrid-scale viscosity. Thus, if the energy balance at the grid scale is written as 

(22) uS 2 {l -Ri t ) =e 

where Ri t denotes an appropriate turbulent Richardson number, and Kolmogorov scaling 

(23) v = Ce 1 ' 3 A 4 / 3 
can be assumed, then Eqs. (22) and (23) give 

(24) v = C 3/2 A 2 Sy/l - Ri t 
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If strong rotation is also present, as in a rotating buoyant flow, Eqs. (17), (18), and (22) imply 

(25) 

Note the difference between the Richardson number dependence in Eq. (24), corresponding to no rota- 
tion, and in Eq. (25), corresponding to strong rotation. Eq. (25) predicts that stabilization by Richardson 
number effects is enhanced by rotation, but that destabilizing Richardson number effects are also enhanced 
by rotation. The second prediction should be tested in numerical simulations. 

Centrifugal instability is another external agency which has often been treated (Bradshaw, 1969) by 
Richardson number modifications of turbulence models. A representative proposal (Launder et al , 1977), 
applicable to flows in which the mean velocity field has circular streamlines 

u r = u z — 0, u$ = U(r) 


is 


(26) Ri t = C R {j) 2 ^(rU) 

where r denotes distance to the rotation axis and + U(r)jr is the total mean angular velocity. In Eq. 
(26), the velocity U should be understood as the velocity relative to the rotation axis: this identification 
maintains the Galilean invariance of the theory. For a subgrid model, it is appropriate to take the turbulent 
frequency scale proportional to the resolved strain rate. With this modification, the turbulent Richardson 
number of Eq. (26) becomes 


(27) Ri t = C' R S~ 2 ^(rU) 

Any part of the mean velocity which corresponds to rigid rotation must be included in Q; equivalently, we 
assume that d 2 U/dr 2 is nonzero. 

Using the formulas for mean vorticity 


W=— - = - — 

dr r r dr 


(rU) 


and mean strain rate 




u 

r 


d.U 

r S < 7> 


we can replace Eq. (27) by the explicitly Galilean-invariant expression 


(28) Rit = \c’ r S- 2 W{W - S) 

The definition of Ri t in Eq. (28) can be substituted in Eq. (25) to provide a subgrid model suitable 
for rotating channel flow or for the problem of rotating Kolmogorov flow recently simulated by Shimomura 
(1995). 

We emphasize that in this treatment of external agencies, Richardson number effects are always production- 
related; they never modify the dissipation rate. Morcoever, stabilization by rotation always corresponds to a 
decrease in the energy transfer by turbulence, hence to a decrease in dissipation rate. A different viewpoint 
is sometimes advanced in engineering modeling: see for example, the survey of models of this type by Rodi 
and Scheuerer (1983). 
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4. The Dissipation Rate Transport Equation. Our earlier work (Rubinstein and Zhou, 1996) 
attempts to implement Leslie’s (1972) suggestion that the dissipation rate transport equation be derived 
from the expression Eq. (10) for energy transfer into the inertial range. A complete treatment would require 
TSDIA in order to evaluate the essentially inhomogeneous diffusion effects. The present account will be 
limited to the production and destruction terms, which are amenable to a homogeneous theory. 

The starting point is then the DIA equations Eqs. (2)-(3) in which the strain rate term 




aUi 

dXm 


— 2k 2 kik p 


dUp 


fiimks 


dU s d 
dx r dk T 


is added to the rotation term. Differentiating the time dependent form of Eq. (10) with respect to time, and 
assuming stationarity of the lowest order TSDIA field, 


^ I ]^ > imn(h){P mrs (p)(5ns( t l)Qir (k) 

Pi rs (k)(2n5 (P)Q rnr (q) + f ds[ 

Jo 

P^rs (p)^m/x(Pi s'jQns (Q) f j ^)Qtr (k, t , s) 

Pjrs (k)Gij (k, s)Qns (p? t ? s)Qmr (<1? f » $) 

"h P^rs (p)^m^(p? ^)^ns(^l) t) ^)Qtr(k, £, s) 

“ Pjra (k, t, s)Q ns (p, £, <s)Q mr (q, £, s) 

+ Pfjrs (P )^m/i (p? f j s'jQns (*1> f j ^)^ir (k, t , s) 

( 29 ) - PjrsMGij (k, t, s)Q ns ( p, t f s)Q mr ( q, <, s)]} 

The time-dependent form of this equation might be of interest in the context of non-equilibrium turbulence 
modeling. Note that the first two terms, distinguished by the absence of any time integration, arise from a 
quasi-normal hypothesis. The remaining terms are corrections due to DIA. 

We follow the program outlined before: to find the dissipation rate transport equation in weakly rotating 
turbulence, we will substitute Kolmogorov scaling forms for the descriptors G and Q; to find this transport 
equation in strongly rotating turbulence, we will substitute descriptors appropriate to strong rotation. 

4.1. The destruction term. The destruction terms D e are those which are independent of strain. 
Consider first the destruction term in non-rotating turbulence. The quasi-normal terms can be shown 
(Rubinstein and Zhou, 1996) to contribute 

(30) = 

when expressed in terms of single-point quantities. The most important conclusion is that the integral 
which leads to Eq. (30) is convergent in the large k limit. This implies that there is no Reynolds number 
dependence in the destruction term. The remaining terms in Eq. (29) lead to the same result. The third 
term, for example, has the form 


[/+ - r]v(pMk,p,q)P(k)P(p)Q(q)Q(k) 

Since the combination tjQ is homogeneous of degree zero, the integral remains convergent, and the form Eq. 
(30) again follows. 

To evaluate the destruction term in strongly rotating turbulence, we observe first that the terms of 
lowest order in Q are the quasi- normal terms which we have already analyzed (Rubinstein and Zhou, 1996). 
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Substituting the energy spectrum with the scaling of Eq. (15) in these terms, there results 

(31) d c = -c%e | n | 

The constant C^ 2 is expressed as a convergent integral. The result of Eq. (31) agrees in the strong rotation 
limit with the rotation correction proposed by Bardina et al (1985). Dimensional analysis obviously cannot 
predict this limit, and other limits for D e have been proposed. The result Eq. (31) is consistent with 
Aupoix’s (1987) finding that the stability of rotating decaying turbulence requires 


D e ~ — | H | a with a < 2 

4.2. The production term. The production terms P e depend on the mean velocity gradient. Begin- 
ning again with turbulence without rotation, we find at once that there can be no production of dissipation 
rate without weak breaking of the isotropy of small scales (Xu and Speziale, 1996). This weak anisotropy 
is introduced using Leslie’s (1972) perturbative DIA theory of shear turbulence. Expanding the single time 
correlation function in a power series in the strain rate, 


Qa (k, t) = (k ,t) + Q t ( 3 1} (k ,t)H — 

where Q (0) is the correlation function of isotropic turbulence and setting G t j = G t -°' the response function 
of isotropic turbulence, Leslie found 


(32) 


Qi]\k,t) = J ds{ G < '°\k,t,s)(--g^-+2kik p k 2 g x ^) < ?r/(^> t > s ) 
+ G^\k > t,s){-^ + 2k j k p k- 2 ^)Q^(k,t,s) 

+ G w {k,t,s)kr Q ( °\k,t,s) 

-Q\f(k,t,s)k r ^ — G (0) (fc,M)} 


Nonzero contributions to production are possible from terms which contain the combination SQ (l) . 
These ter ms will contibute production of dissipation terms proportional to quadratic invariants of the mean 
velocity gradient S^Sij and W tj W l:j where W, 3 is the antisymmetric part of the mean velocity gradient. For 
example, the third term in Eq. (15) contributes 


P c = [J + - /-]SeP(k)P( P )Q ( 1 ) (q)Q ( 0 ) (k) 
(33) = [/ + - J-]5e 2 P(k)P(p)5Q(°>(q)<5( 0 Hk) 


where indices have been dropped. The second power of 0 in Eq. (33) arises from the time integration 
in Eq. (32) required to express in terms of Q^°\ The contribution to P £ of Eq. (33) differs from a 
contribution from the quasi-normal terms only in the factors of 0. Since these factors scale like k~ 2jf3 > they 
do not introduce any high wavenumber divergence, and we find that in terms of single-point quantities, they 
contribute 

(34) Pc ~ KS 2 ~ ^Vt/ 2 

where VU 2 denotes terms quadratic in the mean velocity gradient. These terms must be proportional to the 
invariants SijSij or to W x jW x j 
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The occurrence of terms proportional to W tj W i;j in Eq. (34) would lend theoretical support to the 
procedure of ‘sensitizing to irrotational strains’ introduced by Hanjalic and Launder (1980). Since such 
a term cannot appear in energy production P#, we conclude that Pk and P e may not be related by the 
proportionality 


(35) P e — C e \—P k 

Further evidence against the proportionality of Eq. (35) occurs if the series for P £ is taken to higher 
order in VP by substituting the higher order terms in Leslie’s expansion Eq. (32) in Eq. (29). The result, 

(36) P E = v ~{0{VU) 2 + ^0{SJUf + ■■■} 

can be compared term by term with the result of substituting Yoshizawa’s (1984) expansion of the Reynolds 
stress in the definition of energy production, 


(37) P K = u{0{VU) 2 + — O(VtZ) 3 + • • •} 

There is no reason to anticipate term-by-term equality of the scries in braces in Eqs. (36) and (37). 

Evaluating Eq. (33) using the descriptors Eq. (15) of strongly rotating turbulence leads to the scaling 


(38) P e ~ M-'eVU 2 ~ 2)VC/ 2 

K 

where J'(ft) denotes the rotation-dependent viscosity with the strong rotation limit Eq. (18). 

5. Interpolation Formulas. In principle, the series expansions of i/, D e and P £} in positive powers 
of ft following Shimomura and Yoshizawa (1986) for weak rotation, and in negative powers of ft for strong 
rotation, can be continued to arbitrary order. The problem arises to interpolate rationally between these 
limits to obtain a model valid at intermediate rotation rates. Thus, for the turbulent viscosity, we could 
propose a series of approximations of Pade type for eddy viscosity 


(39) 

v = C v — 

E 

(40) 

r R2 
V = C u 

e 


1 + C l {QK/e) 2S 

1 + D^nK/e ) 2 . 1/2 
1 -f Ci(ftX/£) 2 + C 2 ( fttf /£) 4 ' 


which reduce to the usual eddy viscosity formula with O(ft) corrections for low rotation rates, and to the 
limiting form given by Eq. (18) for strong rotation. The constants C{ and Di could be determined in 
principle by matching to the high and low rotation rate expansions constructed above. The corresponding 
interpolation formulas for D e follow by analogy. 

Interpolation formulas can also be proposed for subgrid scale viscosity including additional production 
or stabilization mechanisms modeled as Richardson number effects. For example, for rotating turbulence 
with added centrifugal stabilization or destabilization, we could set 


(41) v = C s A 2 5Vl -Rit{5 2 + ft 2 /( 1 - Ri t )Y 1/2 

where, repeating the definition given earlier in Eq. (27), 

Rit = C' R S- 2 ~(r u ) - \c' r S~ 2 W{W - S) 
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This formula generalizes the subgrid scale model for strong rotation, Eq. (25), to arbitrary rotation rates. 

Younis (1997) has observed that the lowest order interpolation, Eq. (39), when used with analogs for D e 
and P e does not prove satisfactory in computations. At moderate rotation rates, the reduction in e brought 
about by the modified dissipation rate transport equation causes the leading factor K 2 je in v to increase, 
unless the constant Ci is made sufficiently large. Only when Cl is asymptotically large does the factor in 
braces become small enough to reduce the turbulent viscosity. 

This defect may be due to the property of Eq. (39) that the strong rotation limit fixes the constant C \ , 
which then also determines the rotation correction in the weak rotation limit. It appears that these limits 
are not consistent. The more complex model of Eq. (40) may be more satisfactory, but further investigation 
is essential in order to at least suggest the size of the constants. 

This observation recalls a fact noted earlier: stabilization of turbulence is sometimes treated empirically 
as an increase in the dissipation rate e, because this increase will reduce the (unmodified) turbulent viscosity 
K 2 je. From our viewpoint, the stabilization of turbulence by rotation means that energy transfer is blocked; 
therefore, the dissipation rate e must decrease. But this decrease does not incorrectly enhance the eddy 
viscosity, because the eddy viscosity is also modified by rotation according to Eq. (18). 


6. Conclusions. The extension of the direct interaction approximation to rotating turbulence following 
the suggestions of Zhou (1995) makes possible a rational derivation of a dissipation rate transport equation 
and subgrid scale models for strongly rotating turbulence. Empirical turbulence modeling is unable to 
produce these results because rotation precludes simple dimensional arguments from determining the correct 
form of the modelled equations. 
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